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Abstract 

Two sets of infinitely many exceptional orthogonal polynomials related to the Wilson and Askey- Wilson polynomials 
are presented. They are derived as the eigenfunctions of shape invariant and thus exactly solvable quantum mechanical 
Hamiltonians, which are deformations of those for the Wilson and Askey- Wilson polynomials in terms of a degree £ 
(£ — 1,2,.. .) eigenpolynomial. These polynomials are exceptional in the sense that they start from degree I > 1 and 
thus not constrained by any generalisation of Bochner's theorem. 
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1. Introduction 

Orthogonal polynomials have played important roles 
[11,01 on many different stages of mathematics and physics. 
|For example, the matrix models are one of the recent 
active stages for them in particle physics. As is well- 
known there are infinitely many orthogonal polynomials, 
since they can be constructed by Gram-Schmidt orthonor- 
malisation of the basis {l,x,x , . . .} with respect to any 
given weight (or distribution) dw(x), having finite mo- 
ments J x n dw(x) < oo, n G Z + . A more practical way 
of generating orthogonal polynomials {p n (x)} is through 
the three term recurrence relation (Favard's theorem [3|) 

XPn(x) = A n p n+1 (x)+B n p n (x) + C n p n -i(x) (n > 0), (1) 

with = 0. Here the coefficients A n , B n and C n are 

real and A n —\C n > (n > 1). Conversely all the orthog- 
onal polynomials starting with degree satisfy the above 
three term recurrence relation. Among many orthogonal 
polynomials, a central role is played by those satisfying 
(second order) differential or difference equations, in par- 
ticular those appearing in quantum mechanics or in other 
dynamical problems. According to Bochner all poly- 
nomial solutions of the second order differential equation 



a{x)y" + r(x)y' = X n y, 



(2) 



are classical orthogonal polynomials, that is the Hermite, 
Laguerre, Jacobi and Bessel polynomials. Here, a{x) and 
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t(x) are polynomials in x of deg(cr) < 2, deg(r) < 1. On 
top of these classical orthogonal polynomials, a variety of 
orthogonal polynomials, the so-called (g-)Askey scheme of 
hypergeometric orthogonal polynomials @, H S 0|' are 
known to satisfy second order difference equations with 
pure imaginary or real shifts. They also satisfy restrictions 
due to various generalisation of Bochner's theorem Q. 

In 2008, Gomez-Ullate et al @ introduced two new or- 
thogonal polynomials satisfying second order differential 
equations, which are called exceptional (Ai) Laguerre and 
Jacobi polynomials. They are consistent with Bochner's 
theorem, since the polynomials start with degree 1 rather 
than degree constant term. Later reformulation within 
the framework of quantum mechanics and shape invari- 
ance followed in [it], EH- Recently two sets of infinitely 
many orthogonal polynomials satisfying second order dif- 
ferential equations are introduced by the present authors 
flij ]. These new polynomials could be called Xg Laguerre 
and Jacobi polynomials for each positive integer £. The Xg 
polynomials start with degree £ and lack the members from 
degree to degree £ — 1 but they form a complete basis. 
The £ = 1 examples are those reported earlier jf| [l^, [ll| . 

In this letter we present two new sets of infinitely many 
orthogonal polynomials satisfying second order difference 
equations. They could be called exceptional {Xp) Wilson 
and Askey- Wilson polynomials, since they start with de- 
gree I > 1. They are complete and satisfy various restric- 
tions due to various generalisation of Bochner's theorem 
0]. It is quite natural to expect that these exceptional 
polynomials (Laguerre, Jacobi, Wilson and Askey- Wilson) 
will find fruitful applications in many branches of physics 
and mathematics. The method of deriving these new poly- 



Preprint submitted to Physics Letters B 



November 12, 2009 



normals is, as in the cases of the Xi Laguerre and Jacobi 
polynomials, shape invariance [l3| in quantum mechan- 
ics. These polynomials are obtained as the main part of 
the eigenfunctions of a corresponding Hamiltonian of 'dis- 
crete' quantum mechanics (dQM). Thus the orthogonality 
and completeness of the polynomials are a priori guar- 
anteed. Shape invariance combined with Crum's method 
14[ or the factorisation method 15| or the supersymmet- 
ric quantum mechanics [l6| is a sufficient condition for 
exact solvability of a quantum mechanical Hamiltonian; 
the entire energy spectrum and the corresponding eigen- 
functions can be obtained exactly in an algebraic way. 
Discrete quantum mechanics (dQM) was introduced in 
2004 by the present authors [l7|, [l8| in order to reformu- 
late the theory of hypergeometric orthogonal polynomi- 
als @, S 0] within the framework of quantum mechanics, 
so that various concepts and methods of QM would be 
available for the theory of these orthogonal polynomials. 
This revealed the shape invariance of the Hamiltonians for 
all the hypergeometric orthogonal polynomials [7[, rang- 
ing from the Meixner-Pollaczek, continuous (dual) Hahn, 
Wilson, Askey- Wilson and (g-)Racah together with their 



reductions l£j, 20j . Moreover all these dQM systems turn 
out to be exactly solvable in the Heisenberg picture, too 
211 ] . The positive/negative energy parts of the Heisenberg 
operator solution define the annihilation/creation opera- 
tors, which form dynamical symmetry algebras together 
with the Hamiltonian; for example, the q-oscillator alge- 
bra [22} • Another benefit of dQM is that it provided many 
examples of exactly solvable birth and death processes [23J . 
The Hamiltonians for the two infinite sets of exceptional 
orthogonal polynomials are obtained by deforming those 
for the Wilson and Askey- Wilson polynomials in terms 
of a degree £ eigenpolynomial with twisted arguments, so 
that no singularity would appear in the domains where 
the orthogonality measures are defined. This is essentially 
the same procedure employed for the exceptional Laguerre 



and Jacobi polynomials [12j. The necessary technology, 



the discrete versions of Crum's theorem, was published 
recently Q. 

Here we present our preliminary results on the two sets 
of infinitely many exceptional orthogonal polynomials and 
their discrete quantum mechanical formulation, without 
proof. In the next section, the general setting of dQM 
and the method of shape invariance are recapitulated. In 
section three the relationship between the original shape 
invariance and its multiplicative deformation for the gener- 
ation of exceptional polynomials is explained for two typ- 
ical cases of the Wilson and Askey- Wilson polynomials. 
Sections four and five provide the detailed data of the ex- 
ceptional Wilson and Askey- Wilson polynomials, respec- 
tively. The final section is for a summary and comments. 

2. General setting: shape invariance 

In this letter we consider one dimensional discrete quan- 
tum mechanical systems with pure imaginary shifts. The 



dynamical variables are the real coordinate x (x\ < x < 
X2) and the conjugate momentum p = —id x , which are 
governed by the following generic Hamiltonian: 



H = ^V\x) tf* ^/V*{x) + ^W(x) e" w y/V(x) 

-V(x)-V*(x)=AU, (3) 
A = l t(el p ^/v^xj- e-* p y/V(xj), 
A td = -i(^VWe^ - ^VHxje-l p ). (4) 

It is factorised and positive semi-definite. Here the poten- 
tial function V(x) is an analytic function of x and 7 is a 
real constant to be specified for each specific case, see for 
example (fTTj) - (TT21) for the 7 and the domain x% < x < x-i 
for the two cases discussed in this letter. The ^-operation 
on an analytic function f(x) = ^ n a n x n (a n £ C) is de- 
fined as f*(x) = X)n a n x "' m which a n i s the complex 
conjugation of a n . Since the momentum operator ap- 
pears in exponentiated forms, the Schrodinger equation 
TL4> n (x) — £ n (f> n (x) (n = 0,1,2,...), is an analytic dif- 
ference equation with pure imaginary shifts instead of a 
differential equation. The groundstate is annihilated by 
A, A(f>o = (=> H(f>o = 0) and can be chosen 'real', 
())q(x) = 4>q{x). Throughout this letter we consider those 
systems which have a square-integrable groundstate to- 
gether with an infinite number of discrete energy levels: 
= £ < £1 < £% < ■ ■ ■ ■ 

Shape invariance, a sufficient condition for exact solv- 
ability [l3| , is realised by specific dependence of the poten- 
tial function on a set of parameters A = (Ai, A2, . . .), to be 
denoted by V(x; A), A{\), H(X), £„(A), 4> n {x\ A) etc. The 
shape invariance condition to be discussed in this letter is 

13, m, m, H3i 



A(X)A(\y = kA(\ + 5)U(A + 6) + £x(A), (5) 

in which S is a certain shift of the parameters and k is 
a positive constant. This condition is equivalent to the 
following two relations: 



v(x ~q-\)v*(x 



(G) 



= k 2 V(x; A + 8)V*{x - i T , A + 5), 
V(x + i%;\) + V*(x-i%;\) 

= k(V(x; X + 5) + V*(x;X + 5)) - £i(A). (7) 

Then the entire set of discrete eigenvalues and the corre- 
sponding eigenfunctions of TL = TL{\), 



H(\)4> n (x; A) = £ n (\)<p n (x; A), 



is determined algebraically: 



f n (A) = ^« s S 1 (A + s5), 



(8) 



(9) 



A) oc ^(A)U(A + <$)*••• ^(A + (n - l)5) f 

x o (a;; A + nS). (10) 



2 



3. Deformation of shape invariant systems 



In ii3.ll we recapitulate two known shape invariant sys- 
tems, whose eigenfunctions are described by the Wilson 
polynomial and the Askey-Wilson polynomial. In ^3.21 
we deform them in terms of a degree I eigenpolynomial 
(£ = 1,2,.. .) and present a new shape invariant system. 
The general structure is common for these two cases. Ex- 
plicit expressions of various quantities will be given in fj4] 
and fJSJ respectively. 

3.1. original shape invariant systems 

In [2(| various examples of exactly solvable dQM sys- 
tems are explored with emphasis on shape invariance and 
closure relation. A closure relation is another sufficient 
condition for exact solvability, with which the Heisen- 
berg operator solution for the sinusoidal coordinate r)(x) 
can be obtained exactly [2~ 0| ( see [2t| for further devel- 
opments). The examples in [201 ] are arranged by the sinu- 
soidal coordinates (r](x) = x, x 2 ,cosx) and the Wilson and 
Askey-Wilson polynomials are the most generic members 
for r/(x) = x 2 and cosx, respectively, 

Wilson : i](x) = x 2 , < x < oo, 7=1, (11) 
Askey-Wilson : i](x) — cosx, < x < ir, 7 = logq. (12) 

They contain four parameters A = (Ai, A2, A3, A4) (and one 
more parameter q (0 < q < 1) in the Askey-Wilson case) 
and the data for shape invariance are 



fi 1 1 K — J 

V 2 ' 2 ' 2 ' 2 \ 



1 



Wilson, 
Askey-Wilson. 



(13) 



It is straightforward to verify that the potential function 
V(x;X), (57|) in gland ([71]) in $5j satisfies the shape in- 
variance conditions 

The groundstate is annihilated by A {A<j>o — 0), namely, 



Jv*(x - A) <p (x - i\\ A) 



yjv(x + q;\)ct> (x + q;\), (14) 



and it satisfies also 



<j>o(x; A + 8) = (p(x)^Jv(x + A) (f) (x + i^;X), 
where an auxiliary function ip(x) is defined by 

, s def ( 2x : Wilson, 

I 2 sin a; : Askey-Wilson. 

The eigenfunctions {4> n } have the following form: 

</> n (x;\) = (f> (x;\)P n (r)(x);\), 



(15) 



(16) 



where P n (x;X) is a polynomial of degree n in x and is 
'real', P^ip) = P n {x). In the two cases we are considering, 



P n is the Wilson polynomial or the Askey-Wilson poly- 
nomial, which are orthogonal with respect to the weight 
function (f>o(x; A) 2 : 

4> (x; X) 2 P n (n(x); X) P m (n(x); X)dx = h n (X)S nm . 

(18) 

See g] and fj5] for the explicit forms of the groundstate 
wavefunction 0o(x;A), the definition of the polynomial 
Pn(v( x )'> ^) an d the normalisation constant h n (X). 

The action of the operators A and A* on the eigenfunc- 
tion <p n is 

A(X)4> n (x;X) = f n (X)<(> n -i(x;X + S), (19) 
A(X)Un-i{x; X + 5) = b n ^(X)4> n {x- A), (20) 

in which / n (A) and &„_i(A) are real constants related to 
£n as f„(X)b n -x(X) = S n (X). Their explicit forms are 
given in (f6"T]) in g] and in (I75|) in SJSJ The forward shift 
operator T and the backward shift operator B are defined 
by removing the effects of the groundstate wavefunctions 

JF(A) d = Mx; X + 8)- 1 o A(X) o <£ (x; A) 

= itpix)- 1 (ei p -e~i p ), (21) 

£>(A) d = M^X)- 1 oA(X)i o(j) {x;X + 8) 

= -i(V(x; X)ei p - V*(x; X) e -l p )<p{x), (22) 

and their action on the polynomial P n (r); A) is 

F(X)P n ( V (x);X) =f n (X)P n - 1 (r)(x);X + S), (23) 
B(X)P n ^(ri(x);X + S) =b n -i(X)P n (rj(x);X). (24) 

The difference equation for the polynomials reads 

H(X) = B(X)T(X), (25) 
H(X)P n (ri(x);X) = S n (X)P n (v(x);X), (26) 
or: V(x) (P„ (n(x - 17)) - P n (n{x)) ) 

+ V* (x) (P n (r,(x + * 7 )) - P n (t](x)) ) 

= £ n P n (r)(x)). (27) 

In the last equation (|27p . the explicit parameter A depen- 
dence was suppressed for simplicity. The above formulas 
for the polynomials (|2"Tj) - (f2"T)) are well known [7]]. 

3.2. new shape invariant systems 

By deforming the shape invariant systems in the previ- 
ous subsection in terms of a degree £ eigenpolynomial, we 
present an infinite number of shape invariant systems in- 
dexed by a non-negative integer £ = 0, 1, 2, . . .. For £ = 0, 
it is the original one. Explicit expressions of various quan- 
tities (a^ n) i, a^„,2, ^, n ,ij hi in ) will be given in g] and 



(17) m 



In terms of a twist operation t, 



def 



A= (A 1 ,A 2 ,A 3 ,A 4 ), t(A) = (-A 1 ,-A 2 ,A 3 ,A 4 ), (28) 



we define a degree i polynomial Ct{ x ) from the eigenpoly- 
nomial Pi(x): 



&(x;X) d ^Pt(x;t(X+(£-l)5)). 



(29) 



Then we define a potential function Vi(x; A) and a Hamil- 
tonian Tii(X): 



V e (x;X) = V(x;X + £5) 



Vf(x;\) = V*(x;X + £S) 



£e(r)(x -i^);A) 
, &(v(x - ij); X + 6) 
&(r](xy,\ + 6) 
&(r](x - q);X) 
&(r}(x + q); A) 
&(»j(a: + t7);A + «) 



(30) 



(31) 



&fa(z);A + a) ' 
^(A) ^i(eiPy/v;(x;\) - e~** y/VkfaX)) , (32) 

A(A)t ^-»(V^A)e*" - yfvtfa\)e-**), (33) 

^(A) d = ^(A)U^(A). (34) 

It is interesting to note that the deformation of the poten- 
tial function V(x) is multiplicative in contrast to the addi- 
tive deformation of the prepotential w(x) in the ordinary 
QM [12j . With an appropriate choice of the range of pa- 
rameters, the polynomials & (r)(x); A) and tu(r)(x); A + 6) 
have no zero in the rectangular domain in the complex x 
plane, x\ < Hex < X2, — < Imi < ^ and the Hamil- 
tonian TLi is well-defined and hermitian (self-adjoint). The 
potential function Vt(x) satisfies the shape invariance con- 
ditions ©-©, 



{x-%\)Vl{x-q;X) 



(35) 



= n 2 V e (x; A + S)V e *{x - *7; A + 6), 
Ve(x + q;X) + V;(x-q-,X) 

= K(V e (x; X + 6) + V;{x; A + <5)) - E iA {X). (36) 

By using (fTU]) as a Rodrigues type formula, we obtain the 
complete set of eigenfunctions (n = 0, 1, 2, . . .): 



He(X)4>e, n (x; A) = Ee, n (X)<t>i, n (x; A), 
E e<n (X)=E n (X + £5), 

cf>o(x;X + £8) 



(37) 
(38) 



<l>e,o{x; A) 
ipi(x; A) 



\b(Ti(x-qy,\)\ 

dof 4> {x; X + £5) 



^(rj(x);X + 5), (39) 



\&(r)(x - q);X)\ ' 
<p tt n(x;X) =ipi(x;X)Pi t n(r](x);X). 



(40) 
(41) 



Here a degree ^ + n polynomial in x, the polynomial eigen- 
function Pi >n [x; X) is given by 



dcf 



Pg y n(x;X) = a e , n (x;X)P n (x;X + £5) 

+ b Ln (x-X)P n ^(x;X + £6) 7 (42) 



ai >n (x\ A) = &(af; A + 6) 

+ ai >n<1 (X)&- 1 {x;X + S + S') 

+ a e ^ 2 {X)^ 2 {x;X + 28 + 8'), (43) 

b l>n {x; A) = f 6 w (A)et-i(s; X + S + S'), (44) 

where is 



S' = t(<5) = (—|, —5, 5, 5) 



(45) 



Here the coefficients a^ )rSi x(A), a^,n,2(A) and &^, n ,i(A) are 
real. Although the polynomial eigenfunction P^ in (r](x); A) 
is a degree £ + n polynomial in it has only n zeros in 
the domain x\ < x < x 2 . These polynomials are orthogo- 
nal with respect to the weight function ipi(x; A) 2 : 

2 

ipt(x; X) 2 P e , n (r)(x); X)P^ m (-q(x)\ X)dx = hi t n(X)S nm . 

(46) 

They form a complete basis of the Hilbert space just like 
the Wilson or the Askey- Wilson polynomials P n (rj(x);X) 
in the £ = case. It should be stressed that ipt(x;X) is 
not annihilated by At- 

The action of the operators At and A\ on the eigenfunc- 
tion (j)i >n is 

At(\)4>i,n(x; A) = f t , n {X)4>i in ^i(x; X + 6), (47) 
A t {Xy<t> t , n - 1 {x;X + 8) =bi tn - 1 (\)(j) itn (x;\), (48) 

U,nW = /»(* + = + (49) 

The forward shift operator Tt and the backward shift op- 
erator Bi arc defined in a similar way as before 

T e (X) = ip £ {x; X + 6)- 1 o A e (X) o ^(2; A) 

6(77(0: + ^); A + 5)e*" 



<p(x)&(»7(a;); A) 



-^(^(x-i^A + ^e-*"), (50) 

B/(A) d = Mx; A)" 1 o ^(A)^ o ^(x; X + S) 

= 71 ( 7\+J v(x;\ + e6)Z e (r,(x + q);\)ei'> 
&{r){x);X + 6)\ 

- V*(x; X + tS)tt(v(x - q); \)e-i*)<p(x), (51) 
and their action on the polynomial Pt. n (j]\ A) is 

F t {X)P l>n (ri(xyx) =//,n(A)P/, n -i(r?(a:);A + 5) J (52) 
S^(A)P An _ 1 (r ? (a;);A + 5) = 6^ n _ 1 (A)P An (? 7 (x); A). (53) 

The operator Tii(X) acting on the polynomial eigenfunc- 
tions is defined by 



dcf 



^(A) U =B,(A)^(A) 
= V(x;X + £S) 



x e 



fr(7?(a: + q);A) 

&{v( x - qy a) 

^(r?(x-Z7);A + g) \ 
£,i(v(x); X + S) J 



V*(x;X + U) 



qy,x) 



&(rj(x + i%); A) 

&(»j(s);A + *) 
W/(A)P An (»j(a;);A) = £e, n (A)P?,n ; A) . 



x e 



(54) 
(55) 



4. Exceptional Wilson polynomials 

We present the data in for the Wilson polynomial. 

4-1. original shape invariant system 
We take the four parameters as follows: 



A = f (ax, d2, 03, 04), Re a; > (1 < t < 4), 
{a^, a 2 , ag, 04} = {ai, a 2 , 03, 04} (as a set). 



(56) 



The potential function, ground state, eigenpolynomial, 
etc. are given by 



V(x;X) 



4>o(x; A) 



^ T\ j=1 (a 3 + ix) 
2ix(2ix + 1) : 



dcf 



n 7 =i r ( a i + ix ) 



T{2ix) 



(57) 
(58) 



P n (r7(x);A) = W n (x 2 ; a 1: a 2 , a 3 , a 4 ) = IJ(ai + a 3 -) 
n, n + b — I, a± + ix, a± — ix 



x 4P3 , 

ai + a 2 , ai + 03, a± + a 4 
4 

£ n (A) = n(n + b— l), 6 = a,, 

i=i 

/„(A) = -n(n+&-l), &„(A) = -1, 
27m! (n + 6 - 1)„ 



1), (59) 

(60) 
(61) 



T(2n + b) 



j [ r(n + a J +a fe ). (62) 



i<j<fc<4 



Here W n (x 2 ; 01, aa, 03, a 4 ) is the Wilson polynomial, which 
is symmetric with respect to ai, . . . , a 4 . 

^.D. new shape invariant systems 

We generically restrict the original parameter range (|56p 
as follows: 

01,02 E M, {03, 04} = {03, 04} (as a set), 
0<Oj< Recife (j = l,2;fc = 3,4). (63) 

The polynomial £t(x) is 

A) = W^(a;;-ai-^i,-a2-^ 1 ,a3 + ^i,a4 + ^i), 

(64) 

which is symmetric under ai <-> a 2 and/or 03 «-> 04. It 
is elementary to show that £t(x; A) and &{x\ A + (5) have 
no zero on the half real line < x < 00. The hcrmitic- 
ity of Tii requires a stronger condition of no zero in the 
rectangular domain — ^ < Imi < i. It can be determined 



purely algebraically for each t. For the lowest case, I = 1, 
hermiticity is satisfied by 

(a 3 + a 4 )(aia 2 + 5) < (ai + a 2 )(a 3 a 4 + \). (65) 

The higher £ goes the restriction due to hermiticity be- 
comes less stringent. 

The real coefficients at n ,ii ai >n ^ 2 and bg >n ^i are 



a £,n,l(A) 



ln(a\ + a 2 — &3 — 04 — £ + 1) 



ai + a 2 — 03 — 04 — 2(£ — 1) 

x (ai + a 2 + a 3 + a 4 + 2(n + I - 1)) 1 

x ((ai + f) 2 + (a 2 + f) 2 -(a 3 + f + i - if 

-(a 4 + %+i-l) 2 ), (66) 



(67) 



a*,n,2(A) = JJ(aj - a fe - £ + 1) 

j = l fc=3 

£(£ - l)n(a 3 + a 4 + 2(^-1)) 
(ai + a 2 + n)(oi + a 2 - a 3 - a 4 - 2(£ - 1)) ' 
2 4 

&*,»,i(a) =-nii( a J+ at + ,i +^ 1 ) ( 6s ) 

j'=l fc=3 

£n(a\ + a 2 + n + t — l)( a i +02—03 — 04 — ^+1) 
(ax + a 2 + n)(ax + a 2 + a 3 + a 4 + 2(n + £ — 1)) 

The new normalisation constant h^ n (X) is a simple defor- 
mation of the original one h n (X) |62|) : 



(a 1 +a 2 +n + l)(q 3 + q 4 + n + 21 - 1) 

(ai + a 2 + n)(a 3 + a 4 + n + £ - 1) 
xh n (\ + £S). (69) 



5. Exceptional Askey- Wilson polynomials 

We present the data in f£J] for the Askey- Wilson polyno- 
mial. 

5.1. original shape invariant system 
We take the four parameters as follows: 

q x = f (ai,a 2 ,a 3 ,a 4 ), |a,| < 1 (i = l,...,4), 

{a 1; a 2 , ag, 04} = {ai, a 2 , 03, 04} (as a set), (70) 

where g( A i> A 2, - ) = f (g Al , g^ 2 , . . .). We have one more pa- 
rameter q (0 < 3 < 1) but its dependence is not explicitly 
displayed. The potential function, ground state, eigen- 
polynomial, etc. are given by 



dcf 



(1 - e 2 ")(l - qe 2ix ) 
(e 2 ";q)oo 



</>o(x;A; 

n j= i( a i e ";9)c 

P n (r)(x); A) = p„(cosa;;ai,a 2 ,a3,a4|g) 

= ar"( a l a 2; a l a 3, aiA4; ?)n 



(71) 
(72) 

(73) 



X 403 



'q ™, 010203049" , oie™, oie 
0102, 0103, 0104 



£„(A)-(g-"-l)(l-6o"- 1 ), 6 d ^n«J. (74) 



/ n (A) = gt( g -«_i)(i-6g"- 1 ), 6„(A) 



ft, n (A) = 2tt 



(&g"- 1 ;o)„(&o 2 ";g) c 



(75) 
(76) 



Herep„(x; ai, 02, 03, a^q) is the Askey- Wilson polynomial, 
which is symmetric with respect to a±, . . . , 04. 

5.2. new shape invariant systems 

We generically restrict the original parameter range ([70)) 
as follows: 

ai,02 e M, {03, 04} = {03, 04} (as a set), 
l> aj >\a k \ (i = l,2;k = 3,4). (77) 

The polynomial £f(x) is 

&(a;; A) =pi{x; (axq 2 ) , (a 2 g 2 ) ,a 3 g 2 ,043 2 | j ; 

(78) 

which is symmetric under ai <-> 02 and/or 03 «-> 04. It 
is again elementary to show that £e(x; A) and &(x; A + 6) 
have no zero in the real line segment < x < it. The 
hermiticity of Tie requires a stronger condition of no zero 



in the rectangular domain — ^1 < luix < It can be 
determined purely algebraically for each £ For the lowest 
case, £ = 1, hermiticity is satisfied if 

(ai + a 2 )(l - a 3 a 4 ) - (a 3 + a 4 )(l - aia 2 ) 

> (qi + q-i)( ai a 2 ~ a 3 a 4 ). (79) 



The higher £ goes the restriction due to hermiticity be- 
comes the less stringent. 

The real coefficients a^„ 1, a? >n ^ and b^^ n> i are 

o/,n,i(A) = (« n+ ' _1 (ai + oa) - ^^-^(as + a 4 ) 

+ g £ - 1 (or 1 +a 2 - 1 )-(o 3 : 1 +a 4 - 1 )) (80) 

g?~ l {l - q e )(l - g")(l - a 1 " 1 a 2 ~ 1 a 3 a4^~ 1 )a3a4 

X (1 - oj; 1 a 2 - 1 a 3 a 4 g 2 ^- 1 ))(l - aia 2 a 3 a 4 g 2 (™+^ 1 )) ' 
2 4 

o/,n, 2 (A) = n ri( i °7w _i ) 

j=l fc=3 

(<T £ - 1)(1 - ^(l - q n )(l - a 3 a 4 qW-V) 



(1 - aia 2 g")(l - a^ 1 a 2 l a 3 a 4 q 2 ^-^) 



(81) 



fc/, Bl i(A) = ( 9 - £ - i)a - ? n ) n %-«^ n+ ' _i ) 

j=l fc=3 

(1 - aia 2 (Z n+£ - 1 )(l - ar 1 a^ la 3a4g^ 1 ) 



(1 - ai a 2 q n )(l - a 1 a 2 a 3 a4q 2 ( n+e -V) 



(82) 



The new normalisation constant he in (X) is again a simple 
deformation of the original one h n (X) (f76")) : 



hi,„(X) = q 



_ £ (1 - aia 2 g" +£ )(l - a 3 a 4 <f' 



(1 - aia 2 9 n )(l - a 3 a 4 q n + t - 1 ) 
xh n (X + £S). (83) 

5.5. Zimii to Wilson case 

It is well known that the Wilson polynomials are ob- 
tained from the Askey- Wilson polynomials in a certain 
limit Q together with the corresponding potential func- 
tions, A and A* operators and the Hamiltonians [201 ]. 
This is the reason why various formulas for the excep- 
tional Askey- Wilson polynomials and those of the excep- 
tional Wilson polynomials are common as shown in tj3.2l 
In fact we show that the same limiting procedure will con- 
nect various quantities between the two exceptional poly- 
nomials. Following the Askey- Wilson — » Wilson limiting 
arguments in [20(, let us introduce the rescaled variable 
x w — ^x and set 7 = — j-, q = e~i. This implies 
< x < 7T & < x w < L, p w = fp, e< p = e- pW . 
We write the x-dependence explicitly as TCi — He(x;X), 
Ae = Ag(x;X), etc. By setting A = (af , of , af , af) 
(b = Ylj=i °j V )' the desired limit is achieved by L — ► 00 



L — »oo 



lim (1 - q)- ae &(v(x); A) = H7(v W (x W ); A), (84) 
lim (1 - qr^PtMx); A) = P™ ( V w (x w ); A), 

(85) 

! h+e Mx;X) = ^(x w ;X), (86) 
W *^ w -A), (87) 
(88) 

w,„w. ^ (89) 
(90) 

jf(x W ;A), (91) 

(92) 



lim 

L->oo 

lim (l-q)- 2 V e (x;X) =V e w *(x 

L — >oo 

lim (l-g)-V(x)-^ W (^ W ), 

L — > 00 

lim (l-q)- 1 A e (x;X) = -Af{x w - 1 X), 

L~ >oo 

lim (1 - q)- 2 H e (x; A) = Hf{x w ; A), 

L — >oo 



lim (1 — q)J 7 i(x; A) 

L — > oc 

lim (1 - qy 3 B e {x; A) = -Bf{x w ; A), etc 

L — >oo 

where the superscript w denotes the corresponding quan- 
tity in i j4.2l The orthogonality (|4l)|) also correctly reduces 
to that of UHl 

By construction dQM reduces to the ordinary QM in 
the limit when the shifts become infinitesimal. Thus the 
Jacobi polynomials are obtained from the Askey- Wilson 
polynomials and the Laguerre polynomials from the Wil- 
son polynomials in certain limits [jj, [25[ . It is possible to 
write down the limit formulas connecting various quanti- 
ties of the Xi Askey- Wilson polynomials to those of the 
Jacobi polynomials, as well as those connecting the 
Xg Wilson and Laguerre polynomials. We will defer these 
rather technical results to a later publication. 



6. Summary and Comments 

By deforming two well-known shape invariant discrete 
quantum mechanical systems, corresponding to the Wilson 
and Askey- Wilson polynomials, in terms of their eigen- 
polynomials, two infinite sets of shape invariant systems 
are obtained. Their eigenpolynomials form new types 
of orthogonal polynomials, called Xg Wilson and Askey- 
Wilson polynomials, starting with degree I (I — 1,2,.. .), 
which is the degree of the polynomial used for the de- 
formation. By construction these exceptional orthogonal 
polynomials do not satisfy the three term recurrence rela- 
tion ([T]) but they form complete sets due to the hermiticity 
(self-adjointness) of the Hamiltonian. The allowed param- 
eter ranges are determined by the purely algebraic con- 
ditions that the polynomials used for deformation should 
not have a zero in the rectangular domain x\ < Hex < X2, 
— J^l < Jmx < The explicit ranges for the lowest case 
t = 1 are shown in (p5|) and (|T§|) . The higher I goes the 
restriction due to hermiticity becomes the less stringent. 
It is also possible to find the allowed parameter ranges 
outside of those generically specified ranges (|63|) and (j77|) . 
It would be interesting to try to perform similar defor- 
mation of other shape invariant systems; for example, the 
Meixner-Pollaczek, continuous (dual) Hahn and various re- 
ductions of the Askey- Wilson polynomials. More interest- 
ing would be to formulate the exceptional (q-)Racah poly- 
nomials, the most generic examples of discrete quantum 
mechanics with real shifts. It is a good challenge to clarify 
various properties of these new polynomials, e.g. generat- 
ing functions, the Gram-Schmidt construction, substitutes 
of the three term recurrence relations (fTJ), etc., and to pur- 
sue possible physical applications. 

We thank A. Khare for useful discussion. This work is 
supported in part by Grants-in-Aid for Scientific Research 
from the Ministry of Education, Culture, Sports, Science 
and Technology, No. 19540179. 
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